Search for: All records

Sea-level rise submerges terrestrial permafrost in the Arctic, turning it into subsea permafrost. Subsea permafrost underlies ~ 1.8 million km²of Arctic continental shelf, with thicknesses in places exceeding 700 m. Sea-level variations over glacial-interglacial cycles control subsea permafrost distribution and thickness, yet no permafrost model has accounted for glacial isostatic adjustment (GIA), which deviates local sea level from the global mean due to changes in ice and ocean loading. Here we incorporate GIA into a pan-Arctic model of subsea permafrost over the last 400,000 years. Including GIA significantly reduces present-day subsea permafrost thickness, chiefly because of hydro-isostatic effects as well as deformation related to Northern Hemisphere ice sheets. Additionally, we extend the simulation 1000 years into the future for emissions scenarios outlined in the Intergovernmental Panel on Climate Change’s sixth assessment report. We find that subsea permafrost is preserved under a low emissions scenario but mostly disappears under a high emissions scenario.

 

A key initial step in geophysical imaging is to devise an effective means of mapping the sensitivity of an observation to the model parameters, that is to compute its Fréchet derivatives or sensitivity kernel. In the absence of any simplifying assumptions and when faced with a large number of free parameters, the adjoint method can be an effective and efficient approach to calculating Fréchet derivatives and requires just two numerical simulations. In the Glacial Isostatic Adjustment problem, these consist of a forward simulation driven by changes in ice mass and an adjoint simulation driven by fictitious loads that are applied at the observation sites. The theoretical basis for this approach has seen considerable development over the last decade. Here, we present the final elements needed to image 3-D mantle viscosity using a dataset of palaeo sea-level observations. Developments include the calculation of viscosity Fréchet derivatives (i.e. sensitivity kernels) for relative sea-level observations, a modification to the numerical implementation of the forward and adjoint problem that permits application to 3-D viscosity structure, and a recalibration of initial sea level that ensures the forward simulation honours present-day topography. In the process of addressing these items, we build intuition concerning how absolute sea-level and relative sea-level observations sense Earth’s viscosity structure and the physical processes involved. We discuss examples for potential observations located in the near field (Andenes, Norway), far field (Seychelles), and edge of the forebulge of the Laurentide ice sheet (Barbados). Examination of these kernels: (1) reveals why 1-D estimates of mantle viscosity from far-field relative sea-level observations can be biased; (2) hints at why an appropriate differential relative sea-level observation can provide a better constraint on local mantle viscosity and (3) demonstrates that sea-level observations have non-negligible 3-D sensitivity to deep mantle viscosity structure, which is counter to the intuition gained from 1-D radial viscosity Fréchet derivatives. Finally, we explore the influence of lateral variations in viscosity on relative sea-level observations in the Amundsen Sea Embayment and at Barbados. These predictions are based on a new global 3-D viscosity inference derived from the shear-wave speeds of GLAD-M25 and an inverse calibration scheme that ensures compatibility with certain fundamental geophysical observations. Use of the 3-D viscosity inference leads to: (1) generally greater complexity within the kernel; (2) an increase in sensitivity and presence of shorter length-scale features within lower viscosity regions; (3) a zeroing out of the sensitivity kernel within high-viscosity regions where elastic deformation dominates and (4) shifting of sensitivity at a given depth towards distal regions of weaker viscosity. The tools and intuition built here provide the necessary framework to explore inversions for 3-D mantle viscosity based on palaeo sea-level data.

 

The land surface beneath the Greenland and Antarctic Ice Sheets is isostatically suppressed by the mass of the overlying ice. Accurate computation of the land elevation in the absence of ice is important when considering, for example, regional geodynamics, geomorphology, and ice sheet behaviour. Here, we use contemporary compilations of ice thickness and lithospheric effective elastic thickness to calculate the fully re-equilibrated isostatic response of the solid Earth to the complete removal of the Greenland and Antarctic Ice Sheets. We use an elastic plate flexure model to compute the isostatic response to the unloading of the modern ice sheet loads, and a self-gravitating viscoelastic Earth model to make an adjustment for the remaining isostatic disequilibrium driven by ice mass loss since the Last Glacial Maximum. Feedbacks arising from water loading in areas situated below sea level after ice sheet removal are also taken into account. In addition, we quantify the uncertainties in the total isostatic response associated with a range of elastic and viscoelastic Earth properties. We find that the maximum change in bed elevation following full re-equilibration occurs over the centre of the landmasses and is +783 m in Greenland and +936 m in Antarctica. By contrast, areas around the ice margins experience up to 123 m of lowering due to a combination of sea level rise, peripheral bulge collapse, and water loading. The computed isostatic response fields are openly accessible and have a number of applications for studying regional geodynamics, landscape evolution, cryosphere dynamics, and relative sea level change.

 

The Lofoten and Vesterålen archipelagos are located off the outer coast of northern Norway far from the center of the former Fennoscandian Ice Sheet and near the continental shelf edge. Existing relative sea‐level (RSL) data indicate a pronounced mid‐Holocene transgression and interesting connections with the region's prehistoric human settlement history. Here we present seven new sea‐level index points from isolation basins and five terrestrial limiting points from a coastal sedimentary sequence to refine the region's RSL history. Ingression and isolation contacts in isolation basin sediment cores are identified using sedimentary geochemical data, scanning X‐ray fluorescence profiles and phytoplankton analysis. The ages of these contacts are determined using radiocarbon‐based age models. Our index points range from 11.2 to 1.5k cal a bpand are combined with previously published data to predict the spatiotemporal evolution of sea level in this region using an ensemble of spatiotemporal empirical hierarchical models (STEHME). The new RSL curve constrains the timing of the mid‐Holocene transgression, which occurred from c. 9 to 6k cal a bpwhen sea level increased from −4 to 7 m above present day. From c. 6 to 5k cal a bp, RSL rapidly fell to c. 4 m above present values, and more gradually declined at an average rate of c. 0.8 m ka⁻¹over the last 5k cal a bp.Isobase maps derived using the STEHME show a decrease in the regional shoreline gradients since the transgression maximum from 0.25 to 0.07 m km⁻¹. Our data also better define how RSL variations influenced the location and preservation of coastal settlement locations and harbors from the early Stone Age through historic intervals, improving understanding of regional human–environment interactions.

 

Within extreme continental extension areas, ductile middle crust is exhumed at the surface as metamorphic core complexes. Sophisticated quantitative models of extreme extension predicted upward transport of ductile middle-lower crust through time. Here we develop a general model for metamorphic core complexes formation and demonstrate that they result from the collapse of a mountain belt supported by a thickened crustal root. We show that gravitational body forces generated by topography and crustal root cause an upward flow pattern of the ductile lower-middle crust, facilitated by a detachment surface evolving into low-angle normal fault. This detachment surface acquires large amounts of finite strain, consistent with thick mylonite zones found in metamorphic core complexes. Isostatic rebound exposes the detachment in a domed upwarp, while the final Moho discontinuity across the extended region relaxes to a flat geometry. This work suggests that belts of metamorphic core complexes are a fossil signature of collapsed highlands.

 

It is generally agreed that the Last Interglacial (LIG; ∼130 – 115 ka) was a time when global average temperatures and global mean sea level were higher than they are today. However, the exact timing, magnitude, and spatial pattern of ice melt is much debated. One difficulty in extracting past global mean sea level from local observations is that their elevations need to be corrected for glacial isostatic adjustment (GIA), which requires knowledge of Earth’s internal viscoelastic structure. While this structure is generally assumed to be radially symmetric, evidence from seismology, geodynamics, and mineral physics indicates that large lateral variations in viscosity exist within the mantle. In this study, we construct a new model of Earth’s internal structure by converting shear wave speed into viscosity using parameterisations from mineral physics experiments and geodynamical constraints on Earth’s thermal structure. We use this 3D Earth structure, which includes both variations in lithospheric thickness and lateral variations in viscosity, to calculate the first 3D GIA prediction for LIG sea level. We find that the difference between predictions with and without lateral Earth structure can be meters to 10s of meters in the near field of former ice sheets, and up to a few meters in their far field. We demonstrate how forebulge dynamics and continental levering are affected by laterally varying Earth structure, with a particular focus on those sites with prominent LIG sea level records. Results from four 3D GIA calculations show that accounting for lateral structure can act to increase local sea level by up to ∼1.5 m at the Seychelles and minimally decrease it in Western Australia. We acknowledge that this result is only based on a few simulations, but if robust, this shift brings estimates of global mean sea level from these two sites into closer agreement with each other. We further demonstrate that simulations with a suitable radial viscosity profile can be used to locally approximate the 3D GIA result, but that these radial profiles cannot be found by simply averaging viscosity below the sea level indicator site.